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;  ABSTRACT 

Surface  waves  at  the  water  surface  are  influenced  locally  by 
a  submerged  moving  slender  body.  In  addition  to  the  Kelvin 
wave  pattern  which  would  be  created  by  the  moving  body  when 
travelling  at  uniform  speed  under  smooth  water,  waves  on  the 
surface  can  be  diffracted.;  this  report  is  concerned  with  the 
diffraction  effect.  It  is  assumed  that  the  body  is  moving  at 
constant  speed,  is  not  oscillating,  and  that  the  wave  heights 
are  small  in  comparison  with  the  body  dimensions,  so  that  a 
linearized  theory  can  be  applied.  Using  slender  body  theory, 
relative  wave  heights  can  be  evaluated.  For  monochromatic 
waves,  two  critical,  singular  points,  at  which  the  linearized 
theory  does  not  hold,  occur  when  the  component  of  body  velocity 
in  the  Incident  wave  direction,  and  relative  to  the  wave,  is 
equal  to  the  group  velocity  of  the  wave  with  negative  or  posi¬ 
tive  sign*  at  the  former  point,  disturbances  caused  by  the 
body,  and  existing  ahead  of  the  body*  begin  to  grow  to  infinite 
amplitude;  at  the  latter  point,  monochromatic  waves  behind  the 
body  propagate  at  the  energy  velocity  of  the  wave,  and  are  thus 
also  of  infinite  amplitude.  For  the  coincident  wave-body  direc¬ 
tion  case,  approximate  expressions  for  the  results  are  obtained 
for  three  particular  ranges  of  the  encounter  frequency  and  body 
velocity  parameters.  However,  by  considering  a  wave  spectrum 
with  a  local  spectral  density,  instead  of  a  monochromatic  wave, 
the  singularities  of  the  linear  theory  can  be  removed  and  a  mean 
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increase  in  wave  height  can  he  defined  over  the  entire  encounter- 
velocity  domain.  The  required  integrals  over  the  wave  spectrum 
case  are  discussed;  the  actual  numerical  evaluation  of  these 
integrals  has  ^,et  to  le  carried  out. 

FORMULATION  OF  THE  PROBLEM 

When  a  body  moves  beneath  surface  gravity  waves,  the 
question  arises:  to  what  extent  is  the  propagating  wave  pattern 
affected  by  the  presence  of  the  body?  In  addition  to  the  sta¬ 
tionary  wave  pattern  created  by  the  moving  body  when  travelling 
under  smooth  water,  incident  waves  are  disturbed  by  diffraction; 
this  phenomenon  Is  the  subject  of  the  present  study.  It  is 
assumed  that  the  body  is  slender  and  in  uniform  forward  motion; 
the  initial  discussion  deals  with  monochromatic  surface  waves 
travelling  in  a  fixed  direction. 

The  coordinate  system  is  given  in  Figure  1;  L/2  is  the  unit 
of  length,  L  being  the  length  of  the  body.  The  coordinate  system 
(5»  Rj  C)  is  fixed  in  the  body  and  moves  at  speed  V  in  the  posi¬ 
tive  §  direction.  The  contour  C(C)  of  the  transverse  vertical 
cross  section  is  described  by  the  function  rj  =  f(?,  C).  The 
wave  heights  are  assumed  to  be  small  compared  to  the  body  dimen¬ 
sions  even  when  the  body  is  near  the  free  surface;  consequently, 
a  linearized  theory  can  be  applied. 

Since  the  fluid  can  be  assumed  to  be  incompressible, 
inviscid  and  irrotational,  the  problem  can  be  described  by  a 
velocity  potential  cp,  or,  in  a  dimensionless  notation 
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«PT  =  a  i/gL 


where 


<Pm  =  Y 


cp  + 


co  e 


-lot 


+  cpe 


-iol 


[2] 


with  a  the  frequency  of  wave  encounter,  cpQ  the  steady  part  of 
the  potential  due  to  the  forward  motion  of  the  body,  cp  the 
potential  of  the  incident  wave,  and  cp  characterizing  the  dis¬ 
turbance  of  cp  due  to  the  presence  of  the  body,  cp  ,  cp  ,  and  cp 
w  o  w 

are  all  time  independent,  and  the  encounter  parameter  y  is 
given  by 


v  =  crV/g  [3] 

The  nondimens ional  wave  potential  cp^.  is  given  by 

<p  (S,«|,C)  =  hek(C-2dl  +  lk(5  003  0  +  ”  sin  °>  [4] 

w 

with  a  the  angle  of  the  incident  wave  direction  with  the  negative 
§  axis,  d  the  body  depth  divided  by  L,  and  h  the  wave  height 
parameter  given  by 


h  =  oH/wL  [5] 

with  H  the  wave  height,  and  the  wave  frequency  in  a  co¬ 
ordinate  system  fixed  in  space.  The  non-dimensional  wave 
number  k  is  given  by 
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k  *  u)8L/2g  =  rrL/X 

with  X  the  wave  length.  The  relationship  between  a  and  uj  is 
given  by 


a  =  w  ~  Vt»a  cos  a/g 


By  defining 


6  =  wV/g  =  Pr(2kd)  ,tL  =  V/c  [ 

with  the  Froude  number  based  on  depth 

Fr  *  V/(gLd)1/2  | 

the  wave  encounter  parameter  y  can  be  related  to  the  velocity 
parameter  6  by 


Y  =  8  -  fta  cos  a  [10] 

The  wave  encounter  problem  is  discussed  in  Reference  1 
with  particular  regard  for  the  interpretations  to  be  given  to 
encounter  and  velocity  parameters  so  that  a  wave  potential  of 
the  form  given  by  Equation  [4]  can  be  assumed  to  hold  throughout 
the  o  -  a  plane.  These  results  can  be  summarized  as:  the 
velocity  V  is  always  *0j  for  ft  cos  a  <  1,  the  encounter  fre¬ 
quency  is  equal  to  o  by  Equation  [7],  and  the  encounter  angle 
is  a.;  for  ft  cos  a  >  1,  the  encounter  frequency  is  equal  to  -0, 
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the  encounter  angle  is  a  +  tt,  and  the  wave  encounter  parameter 
given  by  Equation  [lO]  changes  sign  so  as  to  remain  a  0.  These 
interpretations  will  be  discussed  later  when  the  particular 
case  of  coincident  encounter  (a  =  0,  tt)  is  considered. 

There  are  several  solutions  to  the  problem  of  wave  dif¬ 
fraction  over  a  submerged  cylinder,  with  the  waves  travelling 
normally  (Reference  2)  and  obliquely  (Reference  3)  to  an 
infinite  cylinder.  The  solution  given  in  the  present  report 
employs  slender  body  theory  and  assumes  that  the  (finite)  body 
has  a  fine  shape  at  the  bow  and  stern;  it  is  thus  particularly 
well  suited  to  the  case  of  co-incident  wave-body  direction,  but 
not  as  applicable  to  the  normal  encounter  case. 

A  theoretical  treatment  of  the  slender  body  diffraction 
problem  is  given  in  Reference  4,  but  assumptions  restricting 
the  encounter  angle  to  0  and  ir,  and  the  body  velocity  to  very 
small  values,  limit  its  applicability. 

A  detailed  statistical  evaluation  of  surface  wave  diffrac¬ 
tion  is  outlined  in  Reference  5»  but  the  required  scattering 
function  to  which  the  method  would  be  applied  is  not  discussed 
in  any  detail.  The  approach  taken  in  the  present  report  is 
aimed  at,  first,  obtaining  the  scattering  behavior,  and  then 
treating  the  statistics  of  waves  in  a  simplified  manner. 
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THE  GENERAL  SOLUTION 


Following  the  approach  of  the  hydrodynamic  theory  of  ship 
motion,  (References  6,  7)  the  function  cp  may  be  written  as 

i 

<p(I,n,C)  -  J  a?  J“  ac  p(«.c)  o(i  +  (f,,)a  +  (f,  )•)*  [u] 

-1  c  *  ^ 


where  (T»t\pC*)  is  the  point  in  the  ($,h,  C)  co-ordinate  system 
where  the  wave  potential  is  to  be  determined,  F(5,C)  is  the 
unknown  source  strength  at  the  body  surface,  and  G^tlC"?  S,r),C) 
is  the  Green's  function  which  satisfied  the  Laplace  equation, 
the  free  surface  condition,  and  the  radiation  condition.  The 
function  F  is  determined  from  the  boundary  condition  at  the 
body 


an 


[12] 


The  function  G  is  discussed  in  References  6  and  7. 


Approximate  expressions  for  cp  will  be  derived  in  the 
following  paragraphs  for  the  case  of  a  slender  body. 

The  Slender  Body  Approximation 

In  8 lender  body  theory,  the  transverse  dimensions  of  the 
body  are  assumed  to  be  small  with  respect  to  the  length  of  the 
body.  Using  this  assumption,  G  behaves  at  the  body  surface  as 
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3  -  1/((T-S)a  +  (f[-Ti)a  +  (C-C)a)*  +  K(l,s)  [13] 


where  K  Is  a  function  only  of  5  and  §.  Using  Equations  [4]  and 
[ll],  the  boundary  condition  of  Equation  [12]  becomes 

acp(?'c>  •»  +  (C-C)a)*  (1+  (f,c)*)*- 

C 

_  __  -2kd+ik£  cos  a  , 

=  -kh(i  sin  a  sgn(?)  +  (f,-)  sgn  (n))e  /(  1-h( f )a 

[1*] 

Integration  along  the  cross  section  contour  with  respect  to  T 
gives  the  functions 


S(S)  »  J  dC  P(5,C)  (1  +  (f>7)8)* 

C  g 


-  hk  [/3(f)  +  1T(§) 


sin 


■] 


-2kd  +  ik?  cos  a 
e  /2ir 


[15] 


S,|(?)  .  ikS(?)  cos  a  +  T(I)  [l6] 

_  _  -2kd  +  iktf  cos  a 

T($)  »  hk  [P>J  +  sin  aj  e  /2vr  [17] 

where  0(?)  is  the  dimensionless  horizontal  cross-section  offset, 
and  t(5)  the  dimensionless  lateral  vertical  cross-section  offset. 
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The  wave  height  is  determined  from  the  behavior  of  the 
potential  cp  at  the  free  surface.  At  ?  =  2d  the  function  G 
(see  Reference  7)  behaves  as 


G  «  2/((?-?)a  +  r\  3  +  4d! 


>*♦ 


Gj,  +  G 


2  > 


[18] 


with 

-ttG  /2 
n 


An  (q,8)e"2qd+iq^"^  cos  0cos  (tr  sin  9) 
(68qa  cos8  8/k)  +  2yq  cos  8  +  (Yak/6a)-q 


An  (q,8)e“2qd+lq^~^  C0S  0cos  (tr  sin  8) 
(6aqa  cos8  8A)  +  2*yq  cos  6  +  (Yak/68)-q 


B  (q,8)e”2qd+iq^~^  COS  0cos  (tr  sin  9) 

-iL— - - ; -  [19] 

(68q8  cos8  8 A)  -  2Yq  cos  8  +  (Yak/6a)-q 


where 

Ax  -  Bi  -  Y8k/58 

A»  *  -i(  (68q  cos  9A)  +2y) 

B*  ■  i((B*q  cos  ft  A)  -2y)  [20  ] 

and 

■  cos-1  ( 1/4 y)  [21 ] 
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The  integration  contours  M*  and  M»  In  the  q  plane  are  given 
by 


The  poles  ,  qB,  q* ,  and  q4  of  the  Integrands  are 

-  k(l-27  cos  6  +  ( 1-^y  cos  0)*)/268  cos*  0 
q»,4  -  Ic(1+2y  cos  0  +  (1+4y  cos  0)*)/268  cos*  0  f22] 


Assuming  S( 1 )  -  S(-l)  *  0*  the  potential  9  at  the  free 
surface  can  be  written  as 

1 

<p(?»h,2d)  mj d{  2S(«)/((T-4)a  +  t?  +  4d*)^  + 

-1 

1  l 

+/ «  3(«)  Qx  +f  d«  S(?),80,  f*j 

-l  -l 


Since  only  the  wave  pattern  le  of  interest  here*  the  wave 
portion*  9*  of  the  potential  9*  will  be  considered.  The  wave 
contribution  of  0n»  denoted  by  gft*  eaanates  fro®  the  residue 
of  the  four  poles.  Using  aquations  [15]»  [16],  [r]  and  [23]* 
and  Integrating  by  parts*  9  becoaes 
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i  i 

$(§,ti,2d)  »  J  d5  S ( 5) (Si+lk  Ss  cos  a)  +  J  d5  T(?)  g3 


T24] 


where  g^  is  defined,  for  5-5  <  0,  in  front  of  the  body. 


n/ c. 

gn  =  4lJ  d9  An(qx(0),e)e 


-q1(0)(2d-i(5-5)  cos  0) 


X  cos  (inqi  (0)  sin  cos  0)^ 


[25] 


and  for  5-5  >  0,  behind  the  body, 


tt/2  -qa(0)(2d+i(5-T)  cos  0) 

gn  *  4i J  d0  An(qa(e),8)e  cos  (rfla(0)  sin  8)/(1-4y  cos  9)^  + 
0x 

tt/2  -qa(e)(2d-i(5-5)  cos  0) 

+  4i J  dO  Bn(qa(0),0)e  cos  (tr3(8)  sin  8 )/( 1+4-y  cos  0)^- 


t/2  -Q*(0)(2d-i(5-T)  cos  ft) 

-  4lJ"  Bn{q4(0),8)e  cos  (nq*(o)  sin  8)/(1+4y  cos  8)^  [26] 


Prnm  this  it  can  be  concluded  that  one  wave  is  progressing  in 
front  of  the  body  and  the  other  three  are  behind  the  body. 
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The  existence  of  this  wave  pattern  can  easily  be  Illustrated 
by  considering  the  two-dimensional  problem  of  a  line  source  at 
the  surface  of  the  water  (see  References  7  and  8).  It  can  also 
be  understood  by  a  rather  simple  consideration:  In  a  coordinate 
system  fixed  in  space,  gravity  waves  must  satisfy  the  dispersion 
relationship.  Equation  [6], 

u)a  =  2gk/L  or  ca  =  gL/2k  [27] 

In  a  coordinate  system  moving  at  speed  V  in  the  direction  (plus 
or  minus)  of  the  waves  the  frequency  is  given  by 

o  —  2k (c  &  V)/L  *  2gk/L  t  2kV/L  [28] 

This  relation  shows  that  four  waves  exist,  which  have  the  same 
frequency  o,  but  different  wavelengths. 

Range  of  Parameters 

As  an  ' ^lustration  of  the  solution  regimes,  and  the  domains 
of  the  encounter  and  frequency  parameters,  the  case  of  body  and 
wave  in  the  same  direction  will  be  considered.  This  Implies  an 
a  value  of  0  or  ir. 

There  are  then  three  domains  of  interest  covering  the  cases 
of  the  body  heading  into  the  waves,  the  waves  overtaking  the 
body,  and  the  body  overtaking  the  waves.  In  the  first  domain, 
the  angle  a  is  r,  the  body  is  heading  into  the  waves,  6  cos  a 
Is  always  less  than  0,  and 
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Y-B+6*  ,  a-T  [29] 

The  second  domain  occurs  where,  with  a  =  0,  0  s  6  *  1,  and 

Y  -  6  -  8*  ,  a  -  0  [30] 

In  the  third  domain,  with  a  «  0  and  0*1,  the  body  is  over¬ 
taking  the  waves,  and  redefinition  of  encounter  parameters 
implies 

7  -  -5  +  6*  ,  a  -  it  [3l] 

In  the  first  domain,  7  increases  steadily  from  0  as  5 
increases  from  0;  the  point  7  -  1/4,  which  will  become  signif¬ 
icant  in  subsequent  paragraphs,  occurs  at  0  of  ((2)^  -  l)/2. 

In  the  second  domain,  7  increases  from  0  to  a  peak  value  of 
1/4  at  0  *  1/2;  then  decreases  back  to  zero;  It  is  symmetric 
about  0  »  1/2.  The  third  regime  has  7  increasing  with  (5-1) 
in  the  same  way  as  in  region  one. 

APPROXIMATE  SOLUTIONS 

In  order  to  obtain  a  qualitative  solution  for  the  wave  form, 
approximate  evaluations  of  the  integrals  of  Equations  [23]  end 
[26]  can  be  carried  out  where  the  point  of  interest  is  far  from 
the  body  ( |l|  »  1). 

Stationary  Phase 

The  method  of  stationary  phase  may  be  applied  to  the 
Integrals  in  Equations  [25]  and  [26],  as  long  as  all  functions 
except  the  exponential  behave  smoothly  near  the  stationary  point. 
Using  the  series  representations  of  Reference  9#  the  stationary 
point  at  6  •  0  yields  the  results,  for  7  >  1/4,  $-?  <  0, 
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gn>0  [32] 

and  for  l  -  ?  >  ^ 

«n  =  t>  (q3>0,Bn,-1^37r/4)+  ft  (q4,0}Bn,-i,-ir/4)  [33] 


where 

ft  ('im(e),8,Bn,u,v)  = 


ul4(r}V(qm(e),0)e  "m 


-Q  {e)(2d+lu(?-?))+iv 


(2Qm(9)(l-J+U7)^"  |? |) 


¥ 


■[34] 


The  po’es  are  given  by  Equation  [22],  For  y  >  1/4,  qj  and 
qa  have  imaginary  parts  and  thus  produce  an  exponential  decay 
which  results  in  their  omission. 

The  next  term  in  the  g  solution  series  at  the  6=0 

n  -  -  / 

stationary  point  is  of  order  |?j  3/s.  At  the  0  =  ir/2  point, 
however,  a  term  of  order  !?|_1  exists;  Equation  [33]  then  has 
an  additional  term 

i4Bri(q3(7r/2),7r/2)e“2q3^7r,/2^d  ocs  (ryq3  (ir/2)  )/q3  (tr/2)  |5l  [35] 

before  terms  of  l§ l"3/8  appear.  This  result  is  of  Interest 
because  it  shows  that  the  wave  diffraction  phenomenon  is  pre¬ 
dominant  along  the  axis  of  the  body  (r|  =  0)»  with  off-axis 
disturbances  decaying  as  ifj"1.  In  all  subsequent  results, 
however,  only  the  leading  terms  of  order  l?|-1/8  will  be 
carrlf  1. 
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The  wave  diffraction  potential  cp  can  now  be  written,  for 
y  >  1/4  and  ?  >  0,  as 


cp  =  0 


[36] 


while  for  ?  <  0, 


£  =  V(q4(0),-l,-3ir/4)  +  ^(q3  (0), -l,-Mr/4)  [37] 


with 


V(qm(0),u,v)  =  — 


-2d(k+qm(0))  +  iu§qm(0)  +  iv 


(irl5|/2qm(0))^(l-4uY)^ 


-  l(k  cos  a-uq  (0))5 

x  d§(|3+iT  sin  a)e  m 

-1 


[38] 


For  the  case  where  7  <  1/4,  the  poles  qj.  and  q2  contribute, 
leading  to,  for  §-5  <  C, 


gn  =  fi(qi»0,An,  +  1,-tt/4) 


[39] 


and  for  S-5  >  0, 


,+1,+t/4)  +  &(q.»0,Bn, -l,+3ir/4)  + 
+  0(q4 , 0,B  ,  -1,  -rr/4) 


[40] 
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The  wave  diffraction  potential  now  becomes,  for  7  <  1/4  and 


5  >  0, 


V  =  ^(Qi  (0),+1,+tt/4) 


[41] 


and  for  §  <  0, 


5  *  V(qa(0),  +  l,+3ir/4)  +  $(q4  (0), -1-3tt/4)  +  V(q3  (0), -l,+ir/4)  [42] 

The  cross-  ectlonal  area  integral  in  Equation  [38]  for  ^ 
can  be  clarified  by  the  relationships 


d§  0(5)  =  4A/LS 


d?  t(§)  =  4A/La 


[43] 


where  A  is  the  longitudinal-horizontal  cross-sectional  area  of 
the  body,  and  A  the  longitudinal-vertical  cross-sectional  area. 

Several  qualitative  results  can  be  distilled  from  the 
general  results  obtained  so  far.  For  f  ample,  in  the  specific 
case  where  d  is  large  enough  so  that  e~  43  »  e~  48  ,  the  body 

is  moving  directly  Into  the  waves,  with  a  -  tt,  and  1-47  Is  not 
small,  the  term  ^)(qa , -1,  +  f/4)  dominates  the  wave  potential  for 
?  <  0  and  values  of  7  >  1/4.  By  comparing  the  wave  potential 
cp,  at  a  specific  point  ?  =  -2x  behind  the  body,  the  result 


+  TTi/4 


cpw(-2x,2d) 


cpi  (-2x,2d)  *  ae 


[44] 


-tm MwcyoMKW1''  •y-**yal»r*‘  iww'iMwtiw  ewsrp*  iwwwhwhiw  ■  ■  ^  r'"i”'"v-; 
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1s  obtained,  in  which  a  can  be  regarded  as  an  incident  wave 
height  alteration  factor,  given  by 

a  =  4k3/8  e'4dk  A/L8(7r8x8(l+4-y))*  [45] 

This  expression  provides  a  qualitative  picture  of  the 
interference  effect  for  this  particular  case;  the  diffraction 
effect  is  proportional  to  the  inverse  square  root  of  the  dis¬ 
tance  behind  the  body,  is  exponential  with  the  depth  ratio  d, 
proportional  to  the  waterplane  area  and  the  3^2  power  of  the 
dimensionless  wave  number,  and  proportional  to  the  inverse  1/4 
power  of  the  wave  encounter  parameter. 

The  values  of  the  poies  for  the  coincident  wave  problem 
are  shown  in  Figure  2,  and  these  may  be  used  to  assess  the 
relative  values  of  the  contributions  of  the  poles  to  the  total 
solution.  For  the  previous  example,  the  value  of  qsA  is  1 
for  the  entire  range  where  the  body  heads  into  the  waves;  q4A 
is  always  greater  except  when  encounter  y  becomes  large,  where 
its  contribution  approaches  that  of  the  qa  pole.  For  the  over¬ 
taking  range,  the  qs  pole  is  also  smaller  and  so  .it  contributes 
thw  major  portion  of  the  result  excent  y  becomes  large. 

An  expression  similar  to  Equation  [45]  cwuld  be  obtained  for 
this  case  by  placing 

q,(0)A  -  (l+2Y-(l+4Y)*)/2&*-  ((1-0)  /0)*  [46] 


into 
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5  «  ^(qa>  -V+tt/4)  [47] 

and  comparing  the  result  with  Equation  [4]  with  a  =  ir  and 
.5  =  2d. 

A 

The  entire  range  of  values  of  7  has  now  been  covered,  with 
the  exception  of  the  region  where  7  approaches  1/4;  at  these 
points,  the  contributions  from  qj.  and  q2  become  infinite,  but 
In  different  ways  at  different  places.  When  7  >  1/4  but  ap¬ 
proaching,  the  qx  and  qa  contributions,  previously  eliminated, 
must  be  considered  to  obtain  the  way  in  which  these  contribu¬ 
tions  become  singular.  For  7  <  1/4  but  approaching,  the  qx 
and  qa  contributions  become  singular  in  quite  another  way. 

These  two  cases  are  considered  next. 

e  =  47  -1 

The  qx  and  qa  integrals  require  special  consideration  as 
the  stationary  phase  point  approaches  the  integration  interval. 
With  the  small  parameter  c  defined  by 

€  =  47  -  l  [48] 

The  integrals  depending  on  qx  and  qa  may  be  written 

tt/2  _ 

gn  -J  <38  w(9)e1?^9)/(l-47  cos  9)^ 

8 

9a  tt/2  a 

-J  +  /  -J  +  off-*) 

8a  9i 


[*9] 
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where  0!  is  defined  by  Equation  [2l],  For  small  e, 

0i  «  (2e)^  [50] 

and  Equation  [49]  may  be  reduced  to 

(Sa-M*  ifq(0s+9l) 

sn  =  (w(0i)/(y  sin  ex)?) J  de  e  [5l] 

o 

where  a  change  of  variable  has  been  employed.  Evaluation  of 
the  integral  results  in 

gn  «  ±  iBaw(81)e1?qi,a^0x  V(4kYa  f  sin  8X  )  [52] 

which  gives  the  way  in  which  the  g  functions  are  singular  as 
e  -  0. 

The  results  for  the  gn  functions  are,  then,  for  ?-?  <  0 

gn  =  U(qx,0i,An,+l,+Tr/2)  [53] 

where 

ff(qm(e),e,An,u,v)  -  u  Sf  [54] 

with 

5T  -  cos  (tWm(0)  sin  8)6a/4kYa(ir|§|e/qm(0)(l-4Y)^)^  [55] 

and  for  ?-T  >  0 
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gn  =  U  (qa.  9i  >An,+l, -7r/2)  + 

+  ft  (q3,0,Bn,-l,+3ir/4)+  ft  (q4 , 0,Bn» -l,-7r/4)  [56] 

The  wave  potential  function  is  defined 

V(qm(6),uJv)  =  ft(qm(0)iu,v)  X  [57] 

which  is  accurate  for  small  Qj. .  The  wave  potentials  then 
become,  for  ?  <  0, 

*  »  V(qx(0x),+1,T)  [58] 

and  for  ?  >  0, 

5  =  ^(q»(®x  )*+l,0)  +  ft(q4(0),-l,-3ir/4)  +  ft(q*  (0), -1,+t/4)  [59] 
e  =  1  -  4y 

The  qj,  and  qa  contributions  are  also  singular  in  the  four 
regions  where  y  <  1/4  but  approaching  1/4.  These  four  particular 
points  occur  at 

6  -  ((2  -  e)*  -  l)/2  [60] 

for  the  body  heading  into  the  waves,  and  at 

6  »  (1  ±  (*)*)/2 

and 

6  -  ((2  -  e)*  +  l)/2 


[61] 
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for  the  body  overtaking  the  waves 

The  evaluation  of  Equation  [2^]  was  carried  out  numerically 
for  small  values  of  e  to  determine  the  behavior  of  the 
functions.  The  results  are,  for  §-?  <  0, 

gn  =  U(qi,0,An,  +  lJ0)  [62] 

where 

u(qm(e),e,An,u,v)  =  6  X  [63] 

with 

X  =  (4qm(e)e*|?|A)*  m  (ei^T//2/qm(0)  |?l  (irc)^)  [64] 

and  for  S-T  >  0, 

gn  =  U(qa,0,An,  +  l,0)  + 

+  ft(q,,0,Bn,-lf+3Tr/4)+  &(q4,0,Bn> -l,-ir/4)  [65] 

In  a  similar  fashion,  with  the  wave  potential  function  defined 
as 

V(qm(0),u,v)  -  0  X 

the  wave  potential  becomes f  for  T  <  0, 


[66] 
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9  «  V(q1(0),  +  l,+Tr/2)  f67] 

and  for  ?  >  0, 

9  =  V(qa(0),  +  l,+7r/2)  +  ^(q4(0),-l,.3ir/4)  +  $(q3 (0) , -1,  -hr/4)  [68] 

The  representation  of  the  diffraction-altered  wave  field  Is 
now  complete.  The  ordinary  type  of  stationary  phase  representa¬ 
tions  hold  throughout  the  range  of  y  except  there  y  m  1/4,  at 
which  points  either  the  corrections  of  Equations  [54]  and  [57], 
or  of  Equations  [63]  and  [66],  must  be  applied. 

WAVE  SPECTRUM  RESPONSE 

The  results  of  the  previous  sections  can  be  used,  together 

with  an  appropriate  wave -number  spectrum,  to  define  a  mean  wave 

alteration  parameter,  a  ,  In  the  presence  of  a  spectrum  of  waves. 

m 

An  an  example  of  how  this  would  be  carried  out,  the  Incident 
wave  height  reduction  factor  of  Equation  [45]  for  the  head-on- 
encounter  case  will  oe  used.  By  defining  the  parameter 

y  »  2kd  [69] 

Equation  [45]  may  be  rewritten 

a  -  A  y3/2  e"2y/dr.»  {rxd(l+2  Fr  yl/2)/2)1/2  [70] 
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oo 

am  =  A2  b  J  dy  e~4y/(l+2  Fr  y1/2)/irxdL8 

2k  d 
o 

oo 

+  A  b  (2d/irx)1/2J  dy  y“3/2  e"2y/(l+2  Fr  y1/2)1/2  [75] 

2k  d 
o 

Evaluation  of  this  expression  can  be  carried  out  numerically  for 

specific  values  of  x,  d.  A,  k  ,  and  L. 

o 

The  general  problem  which  would,  for  example,  consider  the 
case  of  following  waves,  is  rather  more  difficult;  a  multiple 
term  expression  for  a,  or,  where  y  is  close  to  1/4,  an  expres¬ 
sion  with  logarithmic  components,  is  required.  Since  the  final 

expression  for  a  must  itself  be  evaluated  numerically,  a  general 

m 

numerical  procedure  where  the  appropriate  forms  of  the  potential 
5  are  used  In  the  different  regimes  would  seem  to  offer  the  best 
solution  to  this  complex  problem.  Because  of  the  behavior  of 
the  wave  potentials  near  the  y  »  1/4  singular  points,  and  be¬ 
cause  of  the  logarithmic  and  square  root  nature  of  the  singular¬ 
ities,  integration  over  them  to  obtain  the  result  for  a  wave 
spectrum  can.  In  principle,  be  carried  out;  actual  numerical 
evaluations  are  yet  to  be  done. 

CONCLUSIONS 

The  coincident  wave-body  problem  chosen  to  Illustrate  tha 
diffraction  effect  is  seen  to  possess  a  multiplicity  of  problems 
precluding  a  simple  qualitative  result.  The  case  with  the  body 
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heading  Into  the  waves,  with  y  >  1/4  does,  however,  allow  an 
approximate  result  wh'ch  illustrates  the  way  in  which  a  mean 
wave  height  alteration  can  be  defined. 

Equations  [44]  and  [45]  summarize  these  approximate 
results.  Using  Equation  [25]  for  the  head-on  encounter  case. 
Equation  [45]  can  be  rewritten 

e'8dk/(l+26)  -  aaTrx/l6  (A/L8)8  k3  [76] 

This  relationship  is  shown  graphically  in  Figure  3>  where  values 
of  the  modified  alteration  factor,  the  right  side  of  Equation 
[76],  are  plotted  against  the  quantity  dk  for  values  of  the 
body  velocity  parameter  6  from  0  to  2.5.  As  noted  previously, 
as  6  becomes  large,  the  contribution  from  the  q4  pole  becomes 
significant  for  this  case,  as  shown  in  Figure  2.  Care  must  also 
be  taken,  when  applying  these  results,  to  recall  the  phase  shift 
given  by  Equation  [*’’  ]  when  computing  the  wave  diffraction 
effect,  and  that  the  alteration  factor  a  in  the  same 
equation  multiplies  incident  wave  potential. 
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FIGURE  1  -  COORDINATE  SYSTEM  (Unit  Imgth  .  i  body  Ingil.) 
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FIGURE  2  -  POLE  VALUES  FOR  COINCIDENT  WAVE  PROBLEM 


